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O ' Abstract 



Supersymmetry with bilinear R-parity violation provides a predic- 
tive framework for neutrino masses and mixings in agreement with 
current neutrino oscillation data. The model leads to striking sig- 
. nals at future colliders through the R-parity violating decays of the 

lightest supersymmetric particle. Here we study charged scalar lepton 
decays and demonstrate that if the scalar tau is the LSP (i) it will de- 
cay within the detector, despite the smallness of the neutrino masses, 
(ii) the relative ratio of branching ratios -Br (fi e^Vi) / Br[fi — > 
fJ-J^^i) is predicted from the measured solar neutrino angle, and (iii) 
scalar muon and scalar electron decays will allow to test the consis- 
tency of the model. Thus, bilinear R-parity breaking SUSY will be 
testable at future colliders also in the case where the LSP is not the 
neutralino. 



1 Introduction 



Neutrino physics is one of the most rapidly developing areas of particle 
physics The solar neutrino data, including the recent measurement of 
the neutral current rate for solar neutrinos by the SNO collaboration [0] pro- 
vide strong evidence for neutrino flavour conversion. If interpreted in terms 
of neutrino oscillations, the data indicate a large mixing angle between i/g 
and — Ut-, with a strong preference towards the large mixing angle MSW 
solution (LMA). At 3a one has § 

0.25 < tan^ 00 < 0.83 (1) 

for 1 d.o.f., the best-fit-parameters being 

tan^ 9q = 0.44 Am^ = 6.6 x 10"^ eV^ (2) 

This nicely confirms earlier hints found in ref. 0]. The LMA solution will 
be testable independently by KamLAND and first results are expected 
before the end of the year. In addition, current atmospheric neutrino data 
are most easily explained by z/^ oscillations [^, with the 3a ranges (1 

d.o.f.) 

0.3 < sin^ 9Atm < 0.7 , 1.2 x 10"^ eV^ < Am^ < 4.8 x 10"^ eV^ . (3) 

These data leave little doubt that neutrinos are massive particles after all. 

Unsurprisingly the discoveries in neutrino oscillation physics have trig- 
gered an avalanche of theoretical and phenomenological papers on models 
of neutrino masses and mixings [0, the majority of which are based on one 
variation or the other of the see-saw mechanism p], Here we consider a 

phenomenologically viable alternative, namely, supersymmetry with bilinear 



R-parity breaking terms 1 11 , 12 ] , which in contrast to the seesaw mechanism 
generates neutrino masses at the electro- weak scale. Low-scale schemes for 
neutrino masses have the advantage of being potentially testable in near- 
future accelerator experiments. In this paper we study the implications of 
neutrino physics for charged scalar lepton decays. 

Super symmetric models with explicit bilinear breaking of R-parity (BRpV) 
[ pj] , |12| provide a simple and calculable framework for neutrino masses and 
mixing angles in agreement with the experimental data |13|. BRpV is a hy- 
brid scheme in which one neutrino mass is generated at tree-level, through 
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the mixing with the neutrahnos [|I4[, in an effective "low-scale" variant of 



the seesaw, while the remaining two masses are generated at 1-loop order. 
A complete 1-loop calculation of the neutrino-neutralino mass matrix []T3 



is therefore necessary, before one can confront the model with experimental 
data from atmospheric and solar neutrino experiments. Especially note that 
the "solar" angle has no meaning in BRpV at tree-level. 

BRpV might be considered either as a minimal three-parameter exten- 
sion of the MSSM, valid up to some very high energy scale (like the GUT 
scale) |15| or as the effective description of a more fundamental theory in 
which the breaking of R-parity is spontaneous |T^, 0. While spontaneous 
breaking of R-parity may be considered theoretically more attractive since, 
for example, it provides a motivation for the absence of trilinear R-parity 
breaking parameters in the superpotential, for the sake of simplicity in our 
numerical calculation we will stick to explicit BRpV only. 

One should, however, note that the results obtained here are valid also 
in those classes of models where R-parity is broken spontaneously including 
the presence of an additional Goldstone boson, namely the Majoron J. This 
can be seen as follows: The Majoron consists mainly of the imaginary parts 
of the SU{2) ® U{1) singlet scalars, such as the right-handed sneutrinos [|1^ . 
The only terms which couple the Majoron directly to sleptons is given by 
h'^LH2i^R in the superpotential and the corresponding term in the soft SUSY 
breaking Lagrangian. These terms can in principle induce decays like f 
jl J. However, such a decay requires that one of the charged particles involved 
contains a large left-handed component whereas the other one contains a 
large Higgs component. As we will see below, in the cases we will study the 
sleptons are mainly right-sleptons. In addition, in mSUGRA scenarios the 
mass differences between the lightest three sleptons is rather small leading 
to a further suppression of Majoron-emitting charged slepton decays. 

If R-parity is broken the lightest supersymmetric particle (LSP) will de- 
cay. As was shown in (see also |TB|), if the LSP is the lightest neutralino, 
the measured low-energy neutrino properties translate into predictions for 
the ratios of various branching ratios of the neutralino decay, thus providing 
a definite test of the model as the origin of neutrino masses and mixings. 

However, cosmological and astrophysical constraints on its nature no 
longer apply if the LSP decays. Thus, within R-parity violating SUSY a 
priori any superparticle could be the LSP. In this paper we study the case 
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where a charged scalar lepton, most probably the scalar tau, is the LSP []. 
We calculate the production and decays of f, as well as the decays of e and /i, 
and demonstrate that also for the case of charged sleptons as LSPs neutrino 
physics leads to definite predictions of various decay properties. 

This paper is organized as follows. In Sec. ^we will define the model, dis- 
cuss the charged scalar mass matrix and give some formulas for the two-body 
decays of charged sleptons, which are the most important decay channels. In 
Sec. 1^ we will then discuss production and decays of these particles, with 
special emphasis on possible measurements of R-parity violating parameters. 
Finally, in Sect. ^ we summarize our conclusions. 



2 The model 



Since BRpV SUSY has been discussed in the hterature several times |Tl|, |T2 
p!3| , |2T[| we will repeat only the main features of the model here. We will 
follow the notation of [|1^]. 

The simplest bilinear model (we call it the MSSM) is characterized 
by three additional terms in the superpotential 

W = Wmssm + Wi^p (4) 

where Wmssm is the ordinary superpotential of the MSSM and 

Wijt^ = e^LiHu. (5) 

These bilinear terms, together with the corresponding terms in the soft SUSY 
breaking part of the Lagrangian, 

^soft = ^^o^t^^ + BieiLiHu (6) 

define the minimal model, which we will adopt throughout this paper. The 
appearance of the lepton number violating terms in Eq. (|^) leads in general 
to non-zero vacuum expectation values for the scalar neutrinos (z/j), called Vi 
in the rest of this paper, in addition to the VEVs vu and va of the MSSM 
Higgs fields and H^. Together with the bilinear parameters the Vi 
induce mixing between various particles which in the MSSM are distinguished 



^The case of light stop decays was considered in ref. pC 
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(only) by lepton number (or R-parity). Mixing between the neutrinos and 
the neutrahnos of the MSSM, as mentioned previously, generates a non-zero 
mass for one specific linear superposition of the three neutrino flavour states 
of the model at tree-level. For a complete discussion of 1-loop corrections, 
providing mass for the remaining two neutrino states, see ||T3|| . 

For the decays of the charged sleptons it is necessary to calculate the 
mixings between neutrinos and neutralinos, charginos and charged leptons, 
as well as the charged scalar mixing. Since the various mass matrices can be 
found in [jl3l , we will discuss only the charged scalar mass matrix in the next 
section. 



2.1 The charged scalar mass matrix 

With R-parity broken by the bilinear terms in Eq. (^ the left-handed and 
right-handed charged sleptons mix with the charged Higgs of the MSSM, 
resulting in an (8 x 8) mass matrix for charged scalars. As in the MSSM this 
matrix contains the Goldstone boson, providing the mass of the W-boson 
after electro-weak symmetry breaking. One can rotate away the Goldstone 
mode from this mass matrix, using the following rotation matrix 
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where. 
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(7) 



Wo 



W2 = \Jvl + vl + vl)+V 



(8) 
(9) 
(10) 

(11) 
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This matrix has the property that 



RMl±R^ 



0^ 

ML 

■^7 



(12) 



where M^± is a (7 x 7) matrix and the zeroes in the first row and first column 



5. 

correspond to the (massless) Goldstone state in ^ 
We divide the remaining M^± into two parts, 



gauge. 



ML 



where (Mj-,± 



- m-r + {Mi.r^ (is) 

)^^^] contains only R-parity conserving [R-parity vio- 



lating] terms. Note that in the following we assume for simplicity that there 
is no inter-generational mixing among the charged sleptons. This is moti- 



vated by existing constraints from flavour changing neutral currents and 
is consistent with the minimal SUGRA scenario of the MSSM, which we 
will use in the numerical part of this paper. With this assumption also the 
branching ratio /i — > 67 is small in the bilinear model in agreement with 
experimental data. The R-parity conserving part of ML is given by 

•^7 
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where the dots indicate that the matrix is symmetric and 
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m, 
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i2^'VRp 
9 )^C2/3 
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R3 



(14) 



/2^i?P I 1/7 £\2 2 



(15) 

(16) 
(17) 
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(18) 



with = + ^i)- "^A is the MSSM pseudoscalar Higgs mass parameter 
m\ = {^iB) / {spCp), hf and Ai are the Yukawa couphngs and soft break- 
ing trihnear parameters of the charged lepton of generation z, is the Hig- 
gsino mixing parameter characterizing the superpotential, and C2/3 = cos(2/5), 
where (3 is defined in the usual way as tan (3 = vu/vd- The R-parity violating 
part of M'i± can be written as 



(m: 



2 



'h± {XhlY {^hrY 
Xnn M^(^) 



Am 
X 



RR 



(19) 



The Higgs mass correction and the Higgs-Slepton mixing terms in eq. (p!9|) 
are 



4 2 1 

A 2 '^M2 - 2 "^/S "^i C2/3 2 , n2 2l 

t^D f^_D ■5/3 ^ ^ 



(X^^i)i = —rhl^— - /iei— + ^(/ - (/if )^)t'DWiS/3 



5/3 



(,XHR)i = — ^hfvi{AiSf3 + IJC/3) \=hfeiVD — 



M^^P can be written as, 



(20) 

(21) 
(22) 



M 



2(1) 
LL 



Ami iXLL)l2 {Xll)i3 
(Xii)i2 Aml^ {Xll)23 



LL) 13 



LL) 23 



Am 



L3 



with the diagonal terms given by 

Ami = + + 1(9' + {hfmcl + lig- - g^)^'. 

whereas the off-diagonals are 



(23) 



(24) 



(X^,),2 = ei62 + i-)i-)ml4 



(25) 
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+ ViV2 



/ \ / ^1 \ / ^3 \ 2 2 

{Xll)13^€i€3 + ( — )( — jm^gC^ 



(26) 



+ ViVs 



(^LL)23 = e2e3 + ( )( )"^L,C^ 



(27) 



+ V2Vs 



9")c2f, + ^{hifc2p 



Similarly for mI^^\ 



and 



Am|^ = l(/.f)V-^^'^E^ 



RRJij 



:{hf){hf)v,v, 



(28) 
(29) 



2(1) 

Finally, the matrix Mj^^' has the following peculiar structure, 



2(1) 



LR 



{XLR)n 

(-^L/?)l2 (-^L_r) 22 
(-'^Lr)i3 (-^Li?)23 (-^Lr)33 



where 



LR a 



1 U ' 1) ' r 

(^LR)y = ■7^{hf){ — ){ — )Ci3VD HSfi - AiCf3 

In the above equations we have used the following abbreviation 

_ o o 1 



m- 



m, 



{9' + 9"){vl - V, 



(30) 

(31) 
(32) 

(33) 



With the definitions outlined above, once can easily derive approximate ex- 
pressions for the mixing between the charged Higgs and the charged sleptons 
induced by the R-parity breaking parameters. These are given by 



Note that one expects smOuRi ~ hf smOuii, i-e. the mixing between right- 
handed sleptons and the Higgs should be typically much smaller than the 
left-handed Higgs-slepton mixing. 

Finally, the R-parity conserving mixing between left-handed and right- 
handed sleptons is approximately given by 

sin2^,- =.-,1^^%^. (36) 

" ^Li-^Ri 

2.2 Formulas for two-body decays 

Charged scalar leptons lighter than all other supersymmetric particles will 
decay through R-parity violating couplings. Possible final states are either 
ljUk or qq' . For right-handed charged sleptons {Iri) the former by far dom- 
inates over the hadronic decay mode, since the mixing between 1^ and the 
charged Higgs is small, as explained above. 

In the limit {frtf. , m^,^) <^ one has the simple formula for the two-body 

decays fi /,■ + Uk, 



m? 

p _ Ji 



(OTf f ? + iOZ' f? (37) 



Exact expressions for these couplings can be found, for example, in ref. [12 



Even though in the results presented in this paper we have always calculated 



the couplings appearing in eq. (|37|) exactly using our numerical code, it 
is instructive to consider an approximate diagonalization procedure for the 
various mass matrices. This method is based on the fact that neutrino masses 
are much smaller than all other particle masses in the theory and therefore 
one expects that the bilinear R-parity breaking parameters are (somewhat) 
smaller than the corresponding MSSM parameters. For the charged scalar 
mass matrix all necessary definitions have been given above, for details for 
the corresponding procedure for neutralino and chargino mass matrices we 
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For the case where z 7^ j for l^a Ij I] one finds 
E[(Or?;,^)' + (OS;.^.)T = (-/^fc^^-(^7^^1/i + /^^/,l/2)A,)^(38) 
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+ (^f )^(s/3 sin Ohr, - clsi^Vif 



- (ciMt;)" (39) 

Here cj, = cos^Bj ) and sj, = sm{6j,) where Oj, is the left-right mixing angle 
for /j, sin OnRi characterizes the charged Higgs-(right-handed)-Slepton mixing 
and A is given by 

Aj = eiVD + fiVi. (40) 
The quantities yi and y2 are defined as 



= -^^^ (42) 



2;uDetM 



± 



with DetM^i being the determinant of the MSSM chargino mass matrix. 

While eq. (^) above keeps all R-parity breaking paramters in the expan- 
sion up to second order, eq. ( |5D| ) should be valid in the parameter region in 
which the 1-loop neutrino masses are smaller than the tree-level contribution. 

For the case i = j the corresponding formulas are rather cumbersome and 
therefore of limited utility, except for the case / = e. Here, since /le -C 1 one 
can simplify the couplings to. 



E 

k 



[iOTetJ' + {OZIJ'] ^ 2g''xl\A\' (43) 



The parameter A has been defined above and xi is given by 

X, = ''^'^ (44) 
2DetMxO ^ ^ 

with DetM^o being the determinant of the MSSM neutralino mass matrix 
and M2 the soft SUSY breaking SU (2) mass parameter. 

From eq. (|39|) one expects that various ratios of branching ratios should 
contain rather precise information on ratios of the bilinear R-parity breaking 
parameters, for example, Br{fi eJ2^i)/Br{fi f^J2^i) — {.^i/^2Y- We 
will discuss this important point in more detail in the next section. 
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3 Slepton production and decays 



In this section we will discuss charged slepton production and decay modes. 
In order to reduce the number of parameters, the numerical calculations were 
performed in the mSUGRA version of the MSSM. Unless noted otherwise, 
we have scanned the parameters in the following ranges: M2 from [0,1.2] 
TeV, from [0,2.5] TeV, mo in the range [0,0.5] TeV, Aq/uiq and Bq/uiq 
[-3,3] and tan/? [2.5,10]. All randomly generated points were subsequently 
tested for consistency with the minimization (tadpole) conditions of the Higgs 
potential as well as for phenomenological constraints from supersymmetric 
particle searches. In addition, we selected points in which at least one of the 
charged sleptons was lighter than the lightest neutralino, and thus the LSP. 
This latter cut prefers strongly mo << M2. 

R-parity violating parameters were chosen in such a way [|13| that the 
neutrino masses and mixing angles are approximately consistent with the 
experimental data. A good "fit" to the data would require: a) A^ ~ A,-, 
in order to account for a nearly maximal — z/^ angle in atmospheric 
oscillations, eq. b) Ag < A,-, to fulfil the constraints from z/g-oscillation 
searches at reactors [^; c) |A| ~ [0.05, 2] GeV^, for the atmospheric neutrino 
mass scale, eq. d) ei — 62, to have a large angle in solar oscillations, eq. ^ 
and e) |e|^/|A| ~ [0.1, 10], for the solar mass scale, eq. |. 

In order to investigate the dependence of our results on the assumptions 
about the R-parity violating parameters, we construct three different sets 
of points. Setl was calculated to give an approximate "fit" to the neutrino 
data, as described above. Set2 is similar to Setl, except that ei/e2 has been 
varied in a wider range ([0.1,10]), so as to cover both large and small solar 
angles 0. The last set, called Set3 in the following, is again similar to Setl, 
except that €2/63, which is hardly constrained by neutrino data, is varied in 
the interval £2/63 ~ 0.1 — 2. 

In supersymmetric models in which the scalar leptons have a common soft 
SUSY breaking mass parameter at some high scale (mo in mSUGRA) the 
renormalization group evolution leads to some splitting between the scalar 
taus and the e and jl states at the weak scale. While the lightest mass 
eigenstate in the charged slepton sector is usually mainly the f/j, the eigen- 

^Although at the moment the small angle solar solution is ruled out by a careful analysis 
of the solar data 0], it does not cost us much additional effort to keep this option in mind. 
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100 150 200 250 300 350 400 

mi^ [GeV] 

Figure 1: e^e~ —>■ II production cross section as a function of mi at a Linear 
Collider with 0.8 TeV c.m.s energy. From top to bottom: e (dark, on colour 
printers blue), jl (light shaded, green) and f (dark shaded, red). 



values for cr and jiR are not much heavier, such that also cr and piR decay 
mainly via R-parity violating two-body decays. In our numerical calculation 
we therefore not only consider the decays of tr, but also those of cr and fiR. 
These decays can provide information on the R-parity violating parameters 
not accesible in tr decays and allow for additional cross checks of the con- 
sistency of the model. This is true especially for the case of lepton flavor 
violating slepton decays since from eq. |39| one expects them to be directly 
correlated with the BRpV parameters ej. 

For the calculation of the cross section we have adapted the formulas 



given in p6[ to the bilinear model taking into account correctly all mixing 



effects in the numerical calculation. In Fig. |1| we show the cross section 
(j(e^e~ — > //) in [fb] for -^i = 0.8 TeV as a function of the charged scalar 
mass, for e, ji and f, respectively. Assuming an integrated luminosity of 1000 
fb~^ per year can be achieved at a future linear collider |27, this implies 
that around 10^ scalar muons and scalar taus can be directly produced per 
year. For scalar electrons one expects between 10^ and 10^ produced pairs 
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mi [GeV] 

Figure 2: Charged slepton decay length as a function of mj at a hnear colhder 
with 0.8 TeV c.m.s. energy. From top to bottom: e (dark, on color printers 
blue), p, (light shaded, green) and f (dark shaded, red). 
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per year. Since the three R-parity violating two-body decay channels of the 
right-handed sleptons nearly add up to 100 %, one can expect that individual 
branching ratios will be measured to an accuracy of 1 % if they occur with 
similar strength. 

At the LHC the direct production of right-sleptons is small. As a result, 
they will be produced mainly in cascade decays. The relative cr, jiR and 
yields will depend on the details of the cascade decays involved. Let us con- 
sider for simplicity the case where the cascade decays of the coloured particles 
end up in the lightest neutralino as in the MSSM. Beside the kinematics, the 
resulting number of e/j, jiR and tr arising from these decays depends on the 
nature of the lightest neutralino. When this is mainly bino-like, one expects 
that it decays dominantly into an equal number of e^?, Jxr and f^^'s. As a 
result the number of right-sleptons roughly equal to the number of neutrali- 
nos. Also in case of a wino-like neutralino the amount of cr, jiR and tr will 
be equal. However, in this case the main lightest neutralino decay mode will 
be to a ly-boson and a charged lepton, leaving fewer sleptons to be studied. 
However, as discussed in |1^, in this case the neutralino decay modes can 
be used to probe the large atmospheric neutrino angle. For the case where 
the lightest neutralino is higgsino-like it will decay into a W^-boson and a 
charged lepton, or into a Z-boson and a neutrino, similar to the wino case. 
However for large tan/3 the decay into tr will again be important, even for 
higgsino-like neutralinos. 

In Fig. ^ we show the charged scalar leptons decay length (e, jl and 
f, from top to bottom) as a function of the scalar lepton masses for Set3. 
Very similar results hold for the other sets which are therefore not shown. 
All decay lengths are small compared to typical detector sizes, despite the 
smallness of the neutrino masses. The three generations of sleptons decay 
with quite different decay lengths and thus it should be possible to separate 
the different generations experimentally at a future linear collider. Note that 
the ratio of the decay lengths L(f)/L(/i) is approximately given by {h^/hr)"^. 

As mentioned in the previous section, one expects that ratios of branching 
ratios of various charged slepton decays contain rather precise information 
on ratios of the bilinear parameters e,. That this is indeed the case is shown 
in Fig. ^ for the data of Set2 and in Fig. ^ and Fig. || for the data of Set3. 

As can be seen from these figures, the ratio of charged slepton branching 
ratios are correlated with the ratios of the corresponding BRpV parameters 
ej, following very closely the expectation from Eq. nearly insensitive to 
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Figure 3: Ratios of branching ratios for scalar tau decays (top panel) versus 
{^1/^2)'^, and scalar muon decays (bottom panel) versus (ei/es)^ for Set2. 
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(62/63)^ 

Figure 4: Ratios of branching ratios for scalar electron decays versus (£2/63)^ 
for Sets. 
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variation of the other parameters. Recall, that all the points were generated 
through a rather generous scan over the mSUGRA parameters. Ratios of e^'s 
should therefore be very precisely measurable. Moreover, since only two of 
the three ratios of e^'s are independent it is possible to derive the following 
prediction: 

Br{fi eT,t^i)/Br{fi ^ /iE'^i) : Br{fli eJ2i^i)/Br{fli ^rY^i^i) - 
~ Br{ei ^ iJiT.^i)/Br{ei ^ TY.^i) 
which provides an important cross check of the validity of our bilinear R- 
parity model. Any significant departure from this equality would be a clear 
sign that the bilinear model is incomplete. 

As mentioned in the introduction current solar neutrino data prefer a 
large angle solution (LMA). In the BRpV model the solar angle is mainly 
determined by the ratio ei/e2 [0. A measured solar angle therefore leads to 
a prediction for Br {fi ^ eY^^i)/ Br {fi — > fiJ2^i), as shown in Fig. ^ for the 
data of Set2. With the current limits on tan^ Oq, which are 0.25 < tan^ Oq < 
0.83 for the preferred LMA-MSW solution to the solar neutrino problem 
at 3 0" C.L., one can currently predict that this ratio in the BRpV model must 
be in the range [0.09, 1.8]. Additional input on €2/63, for example 62/63 — 1 
to within 10 % would sharpen the predicted value to [0.15, 1.1]. Obviously, 
also a more precise measurement of the solar angle will lead to a tighter 
prediction in the future. In this context it is worth noting that KamLAND 
[0] should be able to fix the solar angle to within ~ 30 %, if LMA is indeed 
the correct solution to the solar neutrino problem. 

Up to now we have discussed only ratios of R-parity violating parame- 
ters, but charged scalar lepton decays allow, in principle, also to gain infor- 
mation on absolute values of these parameters, as relevant, e. g. to fix the 
scale of neutrino masses determined through the analysis of current solar and 
atmospheric data [^. However, such a measurement would require at least 
some information on MSSM parameters which is at the moment unavailable. 

In Fig. we show the total widths in [eV] for scalar tau decays (top panel) 
and scalar muon decays (bottom panel) for the data of Setl displayed versus 
(|6|//i)^mf. Once and m|- have been measured with some accuracy, one can 
determine the absolute value of |6| from this measurement, provided 62/63 is 
known (for example, from the ratio Br{ei — > fiJ2^i)/Br{ei — >• rj^^i))- 

In a similar way, the decay width of the scalar electron contains informa- 
tion on I A I, as is demonstrated in Fig. ||. A priori knowledge on 62/63 leads 
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Figure 6: Ratios of branching ratios for scalar tau decays versus tan^ ^0 for 
Set2. The top panel shown all data points, the bottom one refers only to 
data points with €2/63 restricted to the range [0.9,1.1]. 
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Figure 7: Total widths in [eV] for scalar tau decays (top) and scalar muon 
decays (bottom) for the data of Setl versus (|e|//i)^mp Once and are 
measured, the widths provide information on the absolute value of |e| = |e|. 
Note that Setl fixes €2/63 ~ 1. In general, this ratio must be known with 
some accuracy, before a value for |e| can be derived from the widths. 



19 



to a tighter correlation, as can be seen from the comparison of the resuhs for 
Set2 and Set3. 

To deduce the value of | A| from this measurement one needs the parameter 
combination defined in Eq. (Q). It contains the MSSM parameters 

Ml, M2, /i and tan/3, which could be determined, for example, if at least 
some of the neutralino and chargino eigenstates are accesible at the LHC or 
a possible linear collider. 



4 Conclusions 

Supersymmetric models with bilinear R-parity breaking provide a simple, 
testable framework for neutrino masses and mixings in agreement with cur- 
rent solar, atmospheric and reactor neutrino oscillation data. The model is 
testable at future colliders if the neutralino is the LSP, as was shown pre- 
viously, as well as in the alternative case where one of the charged scalar 
leptons is the LSP, as we have demonstrated here. 

The measured neutrino mixing angles fix certain ratios of the bilinear 
R-parity breaking parameters and, therefore, lead to well-defined predic- 
tions for the ratio of branching ratios of certain slepton decay modes, which 
should be easily measurable at a future collider such as a high energy lin- 
ear collider. Our main result is shown in Fig. where we display Br {fi 
eJ2^i)/Br{fi fi'J^^i) versus the solar neutrino angle, tan^6'0. 

We have also shown how charged scalar lepton decays allow the determi- 
nation of other parameters of our model, thus providing a definite test that 
bilinear R-parity breaking SUSY is the origin of neutrino masses. 
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